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Honeycomb lattice can support electronic states exhibiting Dirac energy dispersion, with graphene 
as the icon. We propose to derive nontrivial topology by grouping six neighboring sites of honeycomb 
lattice into hexagons and enhancing the inter-hexagon hopping energies over the intra-hexagon ones. 

We reveal that this manipulation opens a gap in the energy dispersion and drives the system into a 
topological state. The nontrivial topology is characterized by the Z 2 index associated with a pseudo 
time-reversal symmetry emerging from the Ce symmetry of the Kekule hopping texture, where the 
angular momentum of orbitals accommodated on the hexagonal ’’artificial atoms” behaves as the 
pseudospin. The size of topological gap is proportional to the hopping-integral difference, which can 
be larger than typical spin-orbit couplings by orders of magnitude and potentially renders topological 
electronic transports available at high temperatures. 

PACS numbers: 03.65.Vf,73.43.-f,73.23.-b 


I. INTRODUCTION 

Honeycomb lattice can host electrons with Dirac-like 
linear dispersion due to its C 3 crystal symmetry^, and 
interests in questing for systems with honeycomb lat¬ 
tice structure flourished since the discovery of graphene 
produced by the scotch-tape technique^r— . The Dirac 
dispersion and the associated chiral property of elec¬ 
tronic wave functions accommodated on honeycomb lat¬ 
tice make it an ideal platform for exploring topological 
states^i^ without external magnetic field. It was shown 
first that a quantum anomalous Hall effect (QAHE) can 
be realized when complex hopping integrals among next- 
nearest-neighboring sites of honeycomb lattice are taken 
into accounifi. Later on it was revealed that the intrinsic 
spin-orbit coupling (SOC) in honeycomb lattice can pro¬ 
vide this complex hopping integrals, which drives spinful 
electrons into a topological state with preserved time- 
reversal (TR) symmetry, known as quantum spin Hall 
effect (QSHE)^“— . Quite a number of activities have 
been devoted towards realizing topological states in elec¬ 
tron systems on honeycomb lattice, such as QAHE by 
strainingi^ii^, twistin g^^d^ and decorating graphenei^, 
and QAHE with spin-polarized edge currents in terms of 
the antiferromagnetic exchange field and staggered elec¬ 
tric potentialii. Honeycomb lattice can also be tuned to 
support topological states in photonic crystal a^^d^ and 
cold atoms^. 

In the present work, we explore possible topological 
properties in honeycomb lattice by introducing a Kekule 
texture in hopping energy between nearest-neighboring 
(NN) sites. We take a hexagonal primitive unit cell 
and view the honeycomb lattice as a triangle lattice of 
hexagons [see the dashed red line in Eig. [iKa)]. When 
the inter-hexagon hopping ti is tuned to be larger than 
the intra-hexagon one to, a topological gap is opened 
at the r point accompanied by a band inversion be- 


(a) © * © © * © 



FIG. 1. (Color online) (a) Honeycomb lattice with a Kekule 
texture in hopping energies between NN sites: to inside 
hexagons as denoted by the green bonds and ti between 
hexagons by red ones. The Red dashed hexagon is the prim¬ 
itive cell of triangular lattice with lattice vectors ai, 02 and 
lattice constant ao = |ai| = |ai|. Numbers 1,...,6 in circle 
index atomic sites within a hexagon, (b) Emergent orbitals 
in the hexagonal artificial atom. 

tween orbitals with opposite spatial parities accommo¬ 
dated on hexagons [see Eig. [11(b)]. A pseudo-TR sym¬ 
metry associated with a pseudospin degree of freedom 
and Kramers doubling in the emergent orbitals are re¬ 
vealed based on Cq point group symmetry, which gener¬ 
ates the Z 2 topology. Eor experimental implementations, 
we discuss that, along with many other possibilities, the 
moleeular graphene realized by placing carbon monoxides 
(CO) periodically on Cu [ 111 ] surface^^ is a very promis¬ 
ing platform to realize the present idea, where the Kekule 
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texture can be controlled by adding extra CO molecules. 

The paper is organized as follows. In Sec. [Ill we start 
with a Cq symmetric Kekule hopping texture among NN 
sites in honeycomb lattice, where a pseudo-TR symmetry 
and pseudospin emerge. In Sec. uni we derive an effec¬ 
tive low-energy Hamiltonian of the system and reveal the 
Z 2 topology. The topological properties of the system is 
demonstrated by the edge state and associated Hall and 
longitudinal conductances in sections [IVl and |Vl We dis¬ 
cuss a promising experimental platform for realizing the 
topological state in Sec. IVIi Concluding remarks are pre¬ 
sented in Sec. ivni 


II. KEKULE HOPPING TEXTURE AND 
EMERGENT ORBITALS 


We start from a spinless tight-binding Hamiltonian on 
honeycomb lattice 

= £0 X 4ci + io y] c\cj + ii X 4'Ci' > (1) 

* (4i) 


with Ci the annihilation operator of electron at atomic 
site i with on-site energy sq satisfying anti-commutation 
relations, (i,jf) and (i', j') running over NN sites inside 
and between hexagonal unit cells with hopping energies 
to and ti respectively [see Fig. Hla)]. The orbitals are 
considered to be the simplest one without any internal 
structure, same as the tt electron of graphene. Below we 
are going to detune the hopping energy ti while keeping 
to constant, and elucidate possible changes in the elec¬ 
tronic state. In this case, the pristine honeycomb lattice 
of individual atomic sites is better to be considered as a 
triangular lattice of hexagons, with the latter character¬ 
ized by Cq symmetry. 

Let us start with the Hamiltonian within a single 
hexagonal unit cell 


Ho = toTt 


/ 0 1 0 0 0 1 \ 
10 10 0 0 
0 10 10 0 
0 0 10 10 
0 0 0 1 0 1 
V 1 0 0 0 1 0 / 


( 2 ) 


where T = [c^, [see Fig. (lla)]. The 

eigenstates of Hamiltonian Hq are given by 

|s) = [l,l,l,l,l,l]^; 

Ip,) = [1,1,0,-1,-1,0]^; 

\Py) = [l,-l,-2,-l,l,2f; 

\d,2_y2) = [l,l,-2,l,l,-2f; 

\d,y) = [1, -1,0, 1,-1, Of 

\fyi3,2_y2^) = [l,-l,l,-l,l,-lf (3) 

with eigenenergies 2to,to,to, —to, —to and —2to respec¬ 
tively, up to normalization factors. As shown in Fig.lHb), 


the emergent orbitals accommodated on the hexagonal 
“artificial atom” take the shapes similar to the conven¬ 
tional d and / atomic orbitals in solids. 

A pseudo-TR symmetry operator can be composed in 
the present system with Cq symmetry: T = UJC with K 
the complex conjugate operator and U = —icTy, where ay 
is the Pauli matrix. It can be checked straightforwardly 
that U corresponds to 7r/2 rotation for p orbitals and 
7r/4 rotation for d orbitals given in Eq. ([3]), which yields 
= —1 in the space formed by the p and d orbitalsi^. 
Therefore, the pseudo-TR symmetry satisfies the relation 
= —1, which is same as that for fermionic particles 
even though the spin degree of freedom of electron has 
not been considered here. This indicates that electrons 
acquire a new pseudospin degree of freedom in the present 
system as far as the low-energy physics is concerned. 

Explicitly the wave functions carrying pseudospins are 
given by the emergent orbitals with eigen angular mo¬ 
mentum 


\P±) = i\Pco) ± i\Pv)) ; \d±) = (| 42 -y 2 ) ± . 


( 4 ) 

Distinguished from the intrinsic spin, the pseudospin 
is directly related to the chiral current density on the 
hexagon. Eor a lattice model, the current density be¬ 
tween two sites is given by Ijk = (i/fi)[tocjc/c — toc\cj]. 
The current distributions evaluated using wave func¬ 
tions in Eq. dl for the pseudospin-up and -down states 
are shown in Eigs. (SJa) and (b) with anticlockwisely 
and clockwisely circulating currents. By considering the 
hexagonal artificial atoms composed by six sites in hon¬ 
eycomb lattice, one harvests states with angular momen¬ 
tum merely from simple orbitals, such as tt electrons in 
graphene. The pseudo-TR symmetry is supported by 
the Co crystal symmetry, sharing the same underlying 
physics with the topological crystalline insulator—. How¬ 
ever, crystal-symmetry-protected topological insulators 
addressed so far need strong SOCs to achieve band inver- 
sions22ri24, which is different from the present approach 
as revealed below. 


III. TOPOLOGICAL PHASE TRANSITION 

We calculate the energy dispersion of Eq. m for sev¬ 
eral typical values of ti (hereafter the on-site energy is 
put as 5o = 0 without losing generality). As shown 
in Eig. [21 there are two two-fold degeneracies at the F 
point corresponding to the two two-dimensional (2D) ir¬ 
reducible representations of Cq point group. Project¬ 
ing the wave functions for ti = 0.9to onto the orbitals 
given in Fig. [I^b), it is found that the topmost two 
valance bands show the character of d orbitals whereas 
the lowest two conduction bands behave like p orbitals 
[see Fig. E^c)], with the order in energy same as those 
listed in Eq. Eor ti = to, the d and p bands become 
degenerate at the F point and double Dirac cones appear 
[see Fig.EJd)], which are equivalent to the ones at K and 
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FIG. 2. (Color online) (a) and (b) Current densities in the 
pseudospin-up channel (p+ or d+) and pseudospin-down chan¬ 
nel {p- or d-) respectively. Band dispersions for the system 
given in Fig. [T] (c) ti = 0.9to (Inset: Brillouin zone of the 
triangular lattice), (d) ti = to and (e) ti = l.lto. Blue and 
yellow are for \p±) and \d±) orbitals respectively, and rainbow 
for hybridization between them. The on-site energy is taken 
sq — 0 . 


K' points in the unfolded Brillouin zone of honeycomb 
lattice with the rhombic unit cell of two sites. When ti 
increases further from to, a band gap reopens at the F 
point. As shown in Fig. He) for ti = l.lto, Ihe valence 
(conduction) bands are now occupied by p {d) orbitals 
around the F point, opposite to the order away from the 
F point, and to that before gap closing. Therefore, a 
band inversion between p and d orbitals takes place at 
the F point when the inter-hexagon hopping energy is 
increased across the topological transition point ti = to, 
namely the pristine honeycomb lattice. 

We can characterize the topological property of the 
gap-opening transition shown in Fig.[2]by a low-energy ef¬ 
fective Hamiltonian around the F point. Since the bands 
near Fermi level are predominated by p and d orbitals, 
it is sufficient to downfold the six-dimensional Hamilto¬ 
nian H(k) associated with the tight-binding model ([T]) 
into the four-dimensional subspace d+,p_, d_]. The 
second term in Eq. m is then simply given by 


/to 0 0 0 \ 

0 -to 0 0 

0 0 to 0 

Vo 0 0 -to J 


( 5 ) 


in the four-dimensional subspace. Contributions from 
the third term in Eq. o to the effective Hamilto¬ 
nian should be evaluated perturbatively. Eirst, we list 
the inter-hexagon hoppings in terms of 6x6 matrices 
ti-i, ti 2 , ti 3 , hi and with 


hi 


/OOOtiOOX 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
Vooo 0 00/ 


/O 0 0 0 0 0 \ 

0 0 0 0 ti 0 
0 0 0 0 0 0 
0 0 0 0 0 0 ’ 
0 0 0 0 0 0 
Voooo 0 0/ 


, tl2 — 


/ooooo 0\ 
0 0 0 0 0 0 
0 0 0 0 0 ti 
0 0 0 0 0 0 
0 0 0 0 0 0 
Vo 0 0 0 0 0 / 


on the basis of [c^^c^]. Eollowing the stan¬ 
dard procedures'^, they can be projected to the subspace 
spanned by d_] 
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With Eourier transformations of matrices in Eqs. 
and m, one obtains the effective low-energy Hamilto¬ 
nian H'{k.) on the basis d+,p_, d_] in the momen¬ 
tum space. Expanding H'(k) around the F point up to 
the lowest-orders of k, one arrives at 


H'{k F) 


H+{k) 0 \ 

0 H.{k) J 


with 


H+{k) = 


H-{k) = 


^aotik^ \ 

— ^aotik- 6t — ^agtik^ / ’ 

^aotik- \ 

-^aotik^ St - ^agtik^ / ’ 


( 7 ) 


( 8 ) 


where St = ti — to, k = {kx, ky), k± = kx ^ iky^ 0 is a 
2x2 zero matrix, and ao is the lattice constant of the 
triangular lattice. For St = 0, the Hamiltonians i^+(k) 
and i7_(k) in Eq. (|8]) are the same as the well-known 
one for honeycomb lattice around the K and K' points^^, 
where the quadratic terms of momentum in the diagonal 
parts can be neglected. 

Eor St > 0, however, the quadratic terms are crucially 
important since they induce a band inversion^!, result¬ 
ing in the orbital hybridization in the band structures de¬ 
noted by the rainbow colors in Eig.[2](e). Associated with 
a skyrmion in the momentum space for the orbital dis¬ 
tributions in the individual pseudospin channels, a topo¬ 
logical state appears characterized by the Z 2 topological 
invarian t It is clear that for St < 0 there is 
no band inversion taking place and thus the band gap is 
trivial as shown in Eig. [2](c). 

It is worthy noticing that, comparing with the Kane- 
Mele model for the honeycomb latticei^iii, the mass term 
St{> 0) in Eq. ([8]) can be considered as an effective SOC 
associated with the pseudospin, namely Aesoc = St. 
Eor St = 0.1 to, a moderate Kakule texture in hopping 
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FIG. 3. (Color online) (a) Band dispersion of a ribbon system 
of 36 hexagons with ti = l.lto cladded from both sides by 10 
hexagons with ti = 0.9to- (b) Real-space distribution of the 
in-gap states associated with the red solid dispersion curves in 
(a), (c) and (d) Real-space distributions of current densities in 
pseudospin up and down channels at the momenta indicated 
by the red and green dots 1 and 2 in (a) within the rhombic 
area sketched by dashed line in (b); the excess currents in 
pseudospin up and down channels are indicated by red and 
green arrows in (b). 


energy, the effective SOC is approximately 3000 times 
larger than the real SOC in magnitude in graphene where 
Asoc — O.lmeV and to = 2.7eV. The huge effective SOC 
is due to its pure electronic character as compared with 
the intrinsic SOC originated from the relativistic effect. 
This is one of the fantastic aspects of the present ap¬ 
proach, which renders a topological gap corresponding 
to temperature of thousands of Kelvin. 


IV. TOPOLOGICAL EDGE STATES AND 
ASSOCIATED CONDUCTANCES 


A. Topological edge states 

In order to check the edge state in the present system, 
we consider a ribbon of hexagonal unit cells of ti = l.lto 
with its two edges cladded by hexagonal unit cells of 
ti = 0.9to. As can be seen in Fig. [31(a), additional states 
appear in the bulk gap as indicated by the red solid curves 
carrying double degeneracy. Plotting the spatial distri¬ 
bution of the corresponding wave functions, we find that 
the in-gap states are localized at the two interfaces be¬ 
tween topological and trivial regions [see Fig. |3l(b)]. As 
displayed in Fig. |3l(c) [(d)], there is an excess upward 
(downward) edge current in the pseudospin-up (-down) 
channel associated with the state indicated by the red 
(green) dot in Fig. |3Ka). 



- 0.1 - 0.05 0 0.05 0.1 

Energy [to] 


FIG. 4. (Golor online) (a) Schematic configuration of a six- 
terminal Hall bar where a topological sample (light blue re¬ 
gion) with ti = l.lto is embedded in a trivial environment 
(gray region) with ti = 0.9to- The size of topological scat¬ 
tering region is 240ao x 120ao, and the width of each semi¬ 
infinite lead is 40ao. The injected current flows along the 
edges of topological sample as indicated by the red parts be¬ 
tween electrodes, (b) Longitudinal and Hall conductances of 
the Hall bar as a function of energy of incident electrons. The 
on-site energy is taken £0 = 0. A rhombic topological sample 
is taken for ease of calculation. 


B. Hall and longitudinal conductances 

At the interface between topological and trivial 
regimes, the crystal symmetry is reduced from Cq to C 2 , 
which breaks the pseudo-TR symmetry in contrast to the 
real TR symmetry. As the result, a mini gap of ^ O.Olto 
[unnoticeable in the scale of Fig.|3l(a)] opens in the other¬ 
wise gapless edge state at the F point due to the coupling 
between two pseudospin channels. In order to quantita¬ 
tively check possible backscatterings caused by this mini 
gap, we perform calculations on the longitudinal and Hall 
conductances based on a Hall bar system as sketched in 
Fig. IHa). It is clear that the current / injected from 
the left electrode divides itself into two branches ac¬ 
cording to the pseudospin states, namely pseudospin up 
(down) electrons can flow only in the upper (lower) edge 
of the Hall bar. By matching wave functions at the in¬ 
terfaces between the six semi-infinite electrodes and the 
topological scattering regio n^^i^^ , one can evaluate the 
transmissions of plane waves scattered among all the six 
leads, and then the longitudinal and Hall conductances, 
Gxx — pxxl{^p\x Pxy) Gxy = Pxyl{,p\x Pxy) 
spectively, by the Landauer-Biittiker formalisut^, where 
Pxx and pxy are the longitudinal and transverse resis¬ 
tances respectively. Similar to the case of QSHE with 
magnetic impurities^, in the present system the values 
of conductivity deviate from the quantized ones when 
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FIG. 5. (Color online) Longitudinal conductance Gxx of the 
topological sample given in Fig. SJa) as a function of the en¬ 
ergy of injected electrons: (a) for several typical values of 
inter-hexagon hopping integrals [1.05to, 0.95to], [l.lto, 0.9to] 
and [1.2to, 0.8to], where the first (second) inside bracket is for 
the topological (trivial) region; (b) with several typical sys¬ 
tem sizes 2Lai x La 2 , where the width of the electrodes is 
fixed at 40ao and the inter-hexagon hopping integral is fixed 
at ti = l.lto and ti = 0.9to for the topological and trivial 
regions respectively. 


Fermi level falls in the mini gap of ^ 0.01 to as shown in 
Fig.EKb). It is noticed, however, that both Gxx and Gxy 
heal quickly after several periods of oscillations that come 
from interferences between the two pseudospin channels. 
It is emphasized that almost perfectly quantized conduc¬ 
tances Gxx = 2e^/ti and Gxy = O^i^ are achieved for 
Fermi level beyond 0.04to up to the bulk gap edge at 
O.lto, where the edge states with almost perfect linear 
dispersions hardly feel the existence of the mini gap and 
essentially no appreciable backscattering exists. 

Now we discuss the hopping integral dependence of the 
longitudinal conductance. The size of scattering region 
is same as in Fig. [41(a) and fixed for all cases. As dis¬ 
played in Fig. [51(a), Gxx saturates to the quantized value 
2e^/h as expected for a Z 2 topological state for all the 
cases with ti = 1.05to, l.lto and 1.2to in the topological 
region (whereas 0.95to,0.9to and 0.8to in the trivial re¬ 
gion correspondingly) when Fermi level is set away from 
the mini gaps, accompanied by oscillations due to inter¬ 
ferences between two pseudospin channels. 

Here we show the sample size dependence of the lon¬ 
gitudinal conductance. We fix inter-hexagon hopping in¬ 
tegrals at l.lto and 0.9to in the topological and trivial 
regions respectively. As displayed in Fig. [51(b), Gxx satu¬ 
rates in all cases to the quantized value 2e^/h when Fermi 
level is shifted away from the mini gap. The topological 



-0.1 -0.05 0 0.05 0.1 

Energyltol 


FIG. 6. (Golor online) (a) Dispersion relations and (b) longi¬ 
tudinal conductances of the topological system same as that 
given in Fig. [dja) except that finite SOG is included. 


edge transport remain unchanged when the topological 
sample becomes large. 


V. REAL SPIN AND QSHE 

In addition to the pseudospin, the real spin degree of 
freedom also contributes to transport properties in re¬ 
alistic systems. In absence of the real SOC, the results 
presented in Fig. [4] remain exactly the same, with an ad¬ 
ditional double degeneracy due to the two real spin and 
thus Gxx = 4e^//i. 

An intrinsic SOC is induced when next-nearest- 
neighbor hoppings in honeycomb lattice are taken into 
accoun t ^^4^ . The low-energy Hamiltonian around the F 
point in Eq. © is then modified to 

m ) ( 9 ) 


with 



H+{k) = ( 

^ —St — vX+ 

6t uX — 

F_(k) = ( 

^ —St + vX+ ^agtik^ 

^ -^aotik+ 

^aohk- 
6t — uX — 


where u = 1 and —1 stand for real spin-up and -down 
states respectively. Therefore, for up real spin SOC en¬ 
hances (suppresses) the topological gap in the pseudospin 
up (down) channel presuming A > 0 [see the left and cen¬ 
tral panels of Fig. Ha)]. As far as A < Jt, the system re¬ 
mains the Z 2 topological state associated with the pseu¬ 
dospin, where electrons with up pseudospin and down 
pseudospin counter propagate at edges, both carrying on 
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FIG. 7. (Color online) Molecular graphene realized by dec¬ 
orating the Cu [111] surface with a triangular lattice of CO 
molecules: (a) ti > to generating topological state, (b) ti < to 
for trivial state. Gray balls are CO molecules decorated by 
STM techniques, and red thick bonds are shorter than green 
thin ones which generates the Kekule hopping textures. 


up and down real spins. The longitudinal conductance 
Gxx saturates to Ae^/h as displayed in Fig. [ni(b). 

When SOC is increased to \ = 5t^ the pseudospin down 
(up) channel with real spin up (down) is driven into a 
semi-metallic state with zero band gap and the Dirac 
dispersion appears at the F point. When SOC goes be¬ 
yond St, this Dirac dispersion opens a gap accompanied 
by a topological phase transition. The system now takes 
a QSHE state where at edges electrons with up real spin 
and pseudospin propagate oppositely to electrons with 
down real spin and pseudospin. Evaluating the longitu¬ 
dinal conductance, one finds that Gxx is quantized ex¬ 
actly to 2e^/h (see Eig. [61(b)), and as shown in the right 
panel of Eig. Ha) there is no mini gap in the edge states, 
as protected by real time-reversal symmetryiSiiii. 


VI. DISCUSSIONS 

Much effort has been devoted towards realizing the 
Dirac-like energy dispersion in artificial honeycomb 
lattices^, ranging from optical lattices^i^ to 2D elec¬ 
tron gases modulated by periodic potential o^^i^^i^^ . All 
these systems provide promising platforms for realizing 
topological properties by detuning effective hopping en¬ 
ergy among NN sites either by modulating muffin-tin po¬ 
tentials or bond lengthes periodically. To be specific, here 
we focus on how to achieve a topological state on the Cu 
[111] surface modulated by triangular gates of carbon 
monoxide (CO) moleculesi^. When extra CO molecules 
are placed at specific positions over the pristine molec¬ 
ular graphene, the bonds of the hexagons surrounding 
them are elongated since the CO clusters enhance local 
repulsive potentials and push electrons away from them, 
which reduces the corresponding electron hopping ener¬ 


gies^. It is extremely interesting from the present point 
of view that Kekule hopping textures have already been 
achieved in experiments^^. We propose to place extra 
CO atoms in the pattern displayed in Eig. [7|(a), where 
the intra-hexagon hopping energy to (green thin bonds) 
surrounding the CO clusters is reduced and the inter¬ 
hexagon hopping energy ti (red thick bonds) is enhanced 
relatively. According to the above discussions, the system 
displayed in Eig. (TKa) with ti > to should take a topo¬ 
logical state. The Kekule hopping texture in Eig. [7|(b), 
dual to that shown in Eig. [7|(a), was realized in exper¬ 
iments^, where the system takes a topologically trivial 
state because ti < to- 

The underlying idea of the present scheme for achiev¬ 
ing the Z 2 topological state is to create artificial orbitals 
carrying on opposite orbital angular momenta and par¬ 
ities with respect to spatial inversion symmetry, and to 
induce a band inversion between them by introducing a 
Kekule hopping texture on honeycomb lattice. In the 
sense that it does not require SOC, the present state 
may be understood as a quantum orbital Hall effect. The 
topological properties can also be extended to photonic 
crystalsi^, cold atoms, and even phonon systems where 
sound waves form band due to periodic configurations 
elastic materials. 


VII. CONCLUSION 

We propose to derive topological properties by modu¬ 
lating electron hopping energy between nearest-neighbor 
sites of honeycomb lattice. Due to the Kekule hopping 
texture with Gq symmetry, atomic-like orbitals emerge, 
which carry a pseudospin degree of freedom characteriz¬ 
ing a pseudo time-reversal symmetry and Kramers dou¬ 
bling. We reveal that the effective spin-orbit coupling 
associated with the pseudospin degree of freedom can be 
larger than the intrinsic one by orders of magnitude be¬ 
cause of the pure electronic origin. The present work of¬ 
fers a new possibility for achieving topological properties 
and related novel quantum properties and functionalities 
at high temperatures. 
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